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Abstract 

Two different types of self-organizing and sustaining ordered motion in fluids or plasmas — one 
is a Benard convection (or streamer) and the other is a zonal flow — have been compared by 
introducing a thermodynamic phenomenological model and evaluating the corresponding entropy 
production rates (EP). These two systems have different topologies in their equivalent circuits: the 
Benard convection is modeled by parallel connection of linear and nonlinear conductances, while 
the zonal flow is modeled by series connection. The "power supply" that drives the systems is also 
a determinant of operating modes. When the energy flux is a control parameter (as in usual plasma 
experiments), the driver is modeled by a constant-current power supply, and when the temperature 
difference between two separate boundaries is controlled (as in usual computational studies), the 
driver is modeled by a constant-voltage power supply. The parallel (series)-connection system tends 
to minimize (maximize) the total EP when a constant-current power supply drives the system. This 
minimum/maximum relation flips when a constant- voltage power supply is connected. 



PACS numbers: 52.25.Fi 



I. INTRODUCTION 



In recent studies on non-equilibrium thermodynamics, the entropy production rate (EP) 
has been highhghted as a promising determinant of self-organized ordered structures. The 
minimum EP (MinEP) principle succeeded to describe various structures self-organized in 
nonlinear dissipative systems [1]. The MinEP principle has its historical origin in Helmholtz' 
minimum dissipation principle that determines the velocity distribution of a stationary vis- 
cous flow [2]. L. Onsager generalized this principle to a variational principle relating thermo- 
dynamic forces (potentials) and fluxes [3, 4]. The MinEP principle proposed by I. Prigogine, 
independently rooted in the study of linear discontinuous systems, has extended the scope 
(by hypothetical studies) toward various nonlinear systems. I. Gyarmati showed the equiv- 
alence of Onsager's variational principle and the MinEP principle in a linear dissipative 
system [5]. A literal application of EP in the variational principle, however, may lead to an 
improper result. For example, consider a heat diffusion on a homogeneous medium where 
Pick's law F = —DVT {F: heat flux, T: temperature, D: diffusion coefficient) applies. 
Then, the total EP is the volume integral J F ■ VT^^dV = J Z)| V logTpdV^, which, how- 
ever, is not the target of minimization leading to Fourier's law (which must be J D\VT\'^dV). 

At the opposite pole to the expanded narrative of MinEP, examples of "maximizing" EP 
have been found in various non-equilibrium systems, suggesting the possibility of antithetical 
"maximum EP (MaxEP)" principle. G. W. Paltridge showed, in his pioneering work [6], 
that the atmospheric heat transfer from the tropical to the polar regions maximizes EP in a 
sense that the temperature contrast tends to be enhanced compared to that of simple heat 
diffusion; see also [7-10]. MaxEP states are often found in turbulent fluid systems. Ozawa et 
al. [1 1] pointed out that fluid-mechanical instabilities play an essential role in increasing EP 
a nonlinear regime. Yoshida and Mahajan [12] formulated a thermodynamic model of heat 
transfer through a self-organizing fluid (plasma), demonstrating the bifurcation of a MaxEP 
operation mode; it is shown that the MaxEP branch is stable when the system is "flux- 
driven," i.e. the heat flux driving the non-equilibrium system is the parameter that controls 
the operation mode; the opposite is a "force-driven" system, in which a thermodynamic force 
(for instance, a temperature difference between two boundaries connected to heat baths) is 
the control parameter. We note that, in usual plasma experiments (as well as in most 
natural processes, like atmospheric heat transfer driven by solar heat source), the input 
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power is controlled (or, given as a determining parameter), which must be transferred to 
some heat sink in a steady state, thus the "flux" is the controlling parameter. On the other 
hand, theoretical analysis or a computational study often uses a force-driven model giving 
a boundary condition on the intensive quantities (like the temperature). Mathematically, 
a fiux-driven condition and a force-driven condition correspond to a Neumann boundary 
condition and a Dirichlet boundary condition, respectively (one may consider a mixed-type 
boundary condition if appropriate). 

Statistical mechanical theories build a foundation for the use of EP as a scalar in vari- 
ational principle: R. C. Dewar [13, 14] related EP to the count of "paths" in phase space 
by invoking Jaynes' formalism [15, 16]. R. K. Niven [17] proposed a different information- 
theoretic evaluation of EP; he also pointed out that the driving condition (what he calls a 
"flow controlled" system is equivalent to the aforementioned fiux-driven system) is essential 
to yield a MaxEP state. 

Here, we note that a "maximum" EP does not mean the maximization of the dissipation 
function in a transient/local process; in fact, transient/local EP is bounded only from below. 
We say maximum or minimum by comparing the total EP of bifurcated operation points 
in (quasi) stationary states. In a general nonlinear system, total EP may be changed by a 
structure (flow in a fiuid/plasma system) on a large scale, even if the ultimate dissipation 
occurs on a small scale (such as the Kolmogorov scale). The MaxEP state is the operating 
mode with the largest EP among possible modes of (quasi) stationary states (other authors 
say "global maximum/minimum of EP" [18, 19]). 

In view of the rich phenomenologies and various theoretical interpretations, it may be 
useful to develop a simple paradigm by which one can compare basic structures as well 
as the modes of operation (or drive) of systems. For this purpose, the problem of heat 
transfer in a plasma provides appropriate examples; a plasma is a strongly nonlinear system 
which exhibits various aspects of non-equilibrium phenomena. As pointed out by Yoshida 
and Mahajan [12], a typical example of bifurcation of a MaxEP state is the H-mode [20], 
which is stably sustained by a given large-enough heat flux. Recently formulated variational 
principle of non-equilibrium thermodynamics [21] can determine the stable operating point 
of a nonlinear system (modeled by a nonlinear impedance), and evaluate either maximum or 
minimum EP depending on how the system is driven (or sustained); the min/max duality 
is described by a Legendre transformation of a generalized dissipation function. 
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In the present study, we discuss not only how the system is driven but also how the non- 
linear impedance is "arranged." In the H-mode, the nonlinearity of the thermal conduction 
is caused by the self-organization of a zonal flow that is directed perpendicular to the heat 
flow; its strong shear suppresses the turbulence, resulting an increased impedance of heat 
transport. However, there is another type of self-organization of an ordered flow — in fluid 
mechanics, known as Benard convection, and in plasma physics, called "streamer" [22] — 
which brings about an "opposite" effect on heat transport; such a flow parallels to the direc- 
tion of heat flow and causes convective heat transport, resulting in a reduced impedance of 
heat transport. While the zonal flow blocks the thermal conduction, the streamer opens up a 
new channel of heat transfer. To model these two different structures by equivalent circuits, 
we consider two different topologies of connecting a nonlinear impedance to a "baseline" 
impedance; one is series connection and the other is parallel connection (here, the baseline 
is determined by the impedance observed when the macroscopic flow is absent, i.e. before 
the self-organization). The blocking (series) impedance changes from zero to a finite value as 
the zonal flow grows and suppresses the turbulent heat transport (the baseline impedance is 
that of the turbulent heat transport). The bypath (parallel) impedance, on the other hand, 
changes from infinity to a finite value as the streamer glows (the baseline impedance is that 
of classical heat diffusion). 

The aim of this study is to describe an abstract (general) balance law that dictates 
bifurcation, stability, and EP of operation points (quasi- stationary macroscopic states) in 
a nonlinear driven system; such a thermodynamic relation will be "mechanism free," hence 
specifying the mechanisms that yield nonlinear impedances or control which type of flow 
self-organizes is not the subject of the present practice. Here we refer to a limited number 
of theoretical arguments on relevant mechanisms of self-organization. The creation of a 
zonal flow is explained by the so-called /3 term in the Rossby wave equation [23] and its 
cousin Hasgawa-Mima equation of drift wave (see a review paper [24]). On the other hand, 
a Benard convection is caused by a buoyancy term measuring a baroclinic effect (which is 
related to an inhomogeneity of entropy caused by the gravity; see [25]). A "bifurcation" of 
the zonal flow and the streamer has been studied for magnetic-curvature driven Rayleigh- 
Taylor instabilities, showing that the strength of dissipation is the control parameter [26] . 
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II. HEAT ENGINE MODEL OF SELF-ORGANIZATION 

We consider a layer bounded on one side by a high temperature core plasma (temperature 
T) and on the other side by a low temperature outer region (heat bath with temperature 
To). The temperature difference (thermodynamic force) and the heat flux are related by a 
thermal impedance Z as T — Tq = ZF. When the system is fiux-driven, a heat flux F is 
given, and then, the inner boundary temperature T is a dependent variable. On the other 
hand, in a temperature-driven system, the temperature difference T — Tq is given, and F is 
a dependent variable. 

We start by reviewing the thermodynamic model of H- mode/zonal flow self- 
organization [12]. In a turbulent heat transport, the thermal impedance Z may be written 
as a sum of two components Z = rjo + rj; rjo represents the background (baseline) turbulent 
heat transport, and rj represents a blocking impedance modeling the suppression of turbu- 
lent heat transport by a zonal flow. The nonlinear impedance 77 is a function of the available 
power P to create the flow. For simplicity, we assume a simple form 

viP) = aP, (1) 

where a is a positive constant modeling the increase of the blocking impedance by creation 
of a zonal flow. The essential nonlinearity of this system is brought about by the power P 
that obeys the thermodynamic second law. In an ideal (quasi- stationary) process (Carnot's 
cycle), the maximum available power a.s F{1 — Tq/T). However, the entropy production in 
the (turbulent) diffusion process diminishes the effective power: denoting = Tq + tiqF, 
we estimate 

The force (temperature difference) and flux are related by the impedance: 

T - To = [r/o + viP)]F. (3) 

The operating point is determined by Eqs. (1), (2) and (3). The equivalent circuit of Eq. (3) 
is shown in Fig. 1 (A). 

Next we formulate a model of streamer (Benard convection). Since a streamer creates 
a new route of heat conduction in parallel to the ambient (baseline) diffusion path, we 
consider a parallel-connection circuit such as Fig. 1 (B). In terms of a baseline conductance 
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FIG. 1: The equivalent circuits of (A) a zonal-flow and (B) a streamer (Benard-convection) models. 

Xo (assumed to be a constant) and a nonlinear conductance x(-P); the flux F and the 
temperature T are related by 

i^=[Xo + x(i^)](T-To). (4) 

The heat flow is (conceptually) divided into two channels; Fl = xo{T — Tq) in the baseline 
conductance and Fnl = x{P){T — Tq) in the nonlinear conductance; see Fig. 1 (B). 
We assume that the nonlinear conductance x{P) is written as 

X{P) = bP + cP\ (5) 

where h and c are constants. The reason why we need the second-order term cP^ will be 
explained later. The power P is determined by (2); here we rewrite it as 

Tq\ r. (. To\ / To 



P = F\^l-^j -Fl(^1-^J = [F-Xo{T-To)]. (6) 

The operating point is determined by Eqs. (4), (5) and (6). 

One may solve the series model Eq. (3) or the parallel model Eq. (4) either for T (with 
respect to a given F) or F (with respect to a given T). Since the series model has been 
analyzed in [12], we focus on the parallel model. Solving Eqs. (4), (5) and (6) for F, we 
obtain two branches: 

Fi={T-To)xo (T<T,) 

^ = ,^^^,3 i^XoT' -{b + AcToXo) ' 
c (T - To) 

+ (1 + 26To + QcT^xo) - [hT^ + ^cT^xo) T + cT^xo] (T > T,) 
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where 

I + 2&T0 + VI + 4&T0 

2b 

is threshold temperature at which the bifurcation occurs. Evidently, F = Fi corresponds to 

baseline linear state, while F = F2 (achieved by P > 0) represents the nonlinear state with 

a streamer; see Fig. 2 (A). 

By these solutions, we find that the critical temperature Tc is determined only by b. The 

bifurcation of the nonlinear solution occurs for every positive b; this is in marked contrast 

to the zonal-fiow (series) model in which the bifurcation needs a sufficiently large a [12]. We 

also find that the nonlinear solution becomes singular if c = 0; a higher order nonlinearity of 

x{P) determines F2. In the temperature-driven system, the nonlinear solution in the regime 

of c > is unphysical, because the condition F2 < Xo{T — Tq) makes the power P negative; 

see Eq. (6). 

The fiux-driven system is more complicated because the nonlinear solution T = T2{F) 
is given as an implicit function. In Fig. 2, we plot the numerical solution of the operating 
points. In the temperature-driven system the nonlinear solution must enhance the heat fiux, 
thus we have to assume c < 0. While the force-driven system has physical solution only for 
c < 0, the fiux-driven system has physical solutions for every c, satisfying T < Tq + F/xo- 
If c < 0, T increases monotonically with F. The temperature difference diminishes when c 
approaches 0, and, T = Tc (constant) for all F when c = 0. With a positive c, T decreases 
monotonically with the increase in F, and approaches Tq. In Fig. 3, we plot the conductance 
Xo + x{P{^jT{F))) of the flux-driven system as a function of F. If c < 0, the nonlinear 
conductance converges to Xo the limit of F — )■ 00. 



III. STABILITY 



We analyze the stability of each branch, and determine which branch is realized. We 
start by the temperature-driven system. We flx the temperature T and deflne the total 
conductance C = Xo + x{P)- The following chain of events is caused by a fluctuation 6F in 
the heat flux: 

dC 

6F-,6C= -6F 

8C 

^ 5F'=iT-To)^6F^a^5F 
7 



(A) (B) 
Force (temperature) driven system Flux driven system 




FIG. 2: The operating points of a nonlinear parallel circuit connected to (A) a constant-force 
(temperature) drive, and (B) a constant-flux drive. The linear solutions are displayed by solid 
lines, and the nonlinear solutions are, by dashed lines. Here we assume Tq = l,xo = = 5 and 
c = (—5,0,10). In (A) the nonlinear branch with c = 10 is unphysical, since the corresponding 
power P becomes negative. 
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FIG. 3: The conductance xo + xiP{FjT{F))) of the flux-driven system, plotted as a function of 
the heat flux F. 

The factor a'^ scales the amplification of fluctuation 6F. The system is stable (unstable) 
if < 1 {a^ > 1). The first step of the events requires certain amount of time while 
other events can occur simultaneously. When At is the time period of one cycle, the flux 
perturbation grows or diminishes as SFn := SF\t=nAt = <y"'SFQ. The amplification factor of 
linear solution is estimated as 

al = b{F - To) (l - ^) • (9) 
On the other hand, the amplification factor of the nonlinear solution is 

. 2 _ (10) 
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FIG. 4: Thermodynamic stability of the flux-driven parallel (streamer) system. The linear branch 
(solid line) is stable in the regime of < -Fc (the critical flux of bifurcation). In i*" > Fc, the 
nonlinear branch (dashed lines) with all c is stable, while the linear branch is unstable. 

Evidently, < 1 for T < T^, while aj > 1 and aj < 1 for T > T^. Therefore, the linear 
solution, which is the unique and stable solution in T < Tc, is destabilized at the bifurcation 
point Tc, and, beyond the bifurcation, the nonlinear branch takes over the linear (smaller 
F) solution as the stable branch. 

Next, we study the stability of the flux-driven system. In this case, it is difficult to 
display the amplification factor explicitly. Here we consider the chain of events caused by a 
fluctuation 6T in the temperature: 



We plot the numerical estimate of the amplification factor of each branch in Fig. 4. The 
linear solution destabilizes beyond the critical flux Fc, and the bifurcated nonlinear branch 
is stable for every value of c. 

IV. COMPARISON OF EP IN PARALLEL AND SERIES SYSTEMS 

Here, we compare EP in parallel/series systems connected to a flux/temperature drives. 
The entropy production caused by heat flow is given by the volume integral J -F- V(l /T)dV, 
where F is the local heat flux. In a stationary system, V ■ .F = 0, thus, integrating by parts, 



6T 6Z 




-> 6T' 




9 



we may represent EP by the boundary terms as F{\/Tq — 1/T). Figure 5 shows EP in 
the flux-driven and temperature-driven systems. In the series-connection system, EP is 
maximized by the nonhnear solution when it is flux-driven (A), but is minimized when it is 
temperature-driven (B). On the contrary, EP of the parallel-connection system is minimized 
by the nonlinear solution when it is flux-driven (C), and is maximized when it is temperature- 
driven (D). R. K. Niven [18] studied numerically the EP of turbulent flow in parallel pipes 
driven by a fixed head (force), and found similar behavior. 




— linear branch 
nonlinear branch 



FIG. 5: The entropy production rate (EP) of the series/parallel-connection systems in flux- 
driven/temperature-driven operations, which is either maximized or minimized, by the nonlinear 
branch of solution (that is always the stable branch when it exists), depending on the topology of 
connection (series or parallel), as well as how the system is driven. 

Table I summarizes the EP of the series-connection and the parallel-connection systems 
in flux-driven and temperature-driven operations. 



V. CONCLUSION 



We have introduced a paradigm of describing operating points (macroscopic quasi- 
stationary state) of a nonlinear heat conduction system; two different connections (series 
and parallel connections) of nonlinear impedances (modeling the self-organization of zonal 
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TABLE I: EP of the series-connection and the parallel-connection systems in flux-driven and 
temperature-driven modes of operation. 





Flux-driven 


Temperature-driven 


Series-connection system 


Max 


Min 


Parallel-connection system 


Min 


Max 



flow and streamer, respectively), as well as two different types of drives (flux drive and 
temperature drive), frame the paradigm. Thermodynamic stability of bifurcated operating 
points can be understood by the phenomenological model based on this paradigm; the non- 
linear branch (created by a positive power P) is always stable (when it bifurcates from the 
linear branch). The parallel and series-connection systems, respectively, minimize and max- 
imize EP when connected to a flux drive. When they are connected to a temperature-drive, 
however, the minimum/maximum relation flips. 
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